Domain wall heat conductance in ferromagnetic wires 
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We present a theoretical study of heat transport in electrically insulating ferromagnetic wires 
containing a domain wall. In the regime of validity of continuum micromagnetism a domain wall 
is found to have no effect on the heat conductance. However, spin waves are found to be reflected 
by domain walls with widths of a few lattice spacings, which is associated with emergence of an 
additional spin wave bound state. The resulting domain wall heat conductance should be significant 
for thin films of Yttrium Iron Garnet with sharply defined magnetic domains. 
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Spin wave (SW) excitations in magnetic systems have 
been studied for many decades, but due to recent progress 
in the fabrication of magnetic nanostructures, novel de- 
tection techniques, and new discoveries such as current- 
induced magnetization dynamics, the field is very much 
alive It has been shown that SW can deliver signal 
information and probe magnetic properties. Several de- 
vice concepts have been proposed to use SW logic devices 
d, Q . Recent experiments show that three crucial mate- 
rial parameters, namely, the spin polarization P, the in- 
trinsic Gilbert damping a, and the dissipative correction 
to adiabaticity /3, can be simultaneously determined by 
measuring the current-induced SW Doppler effect 0, [a ■ 
Non- volatile data storage devices @ and logic circuits [3| 
make use of magnetic domain walls (DWs). The combi- 
nation of both strategies is promising as well 0]. Sev- 
eral phenomenon based on the interaction between SWs 
and DWs have been proposed within the framework of 
continuum micromagnetic theory, such absence of reflec- 
tion [§] but an associated scattering phases shift [1, Q 
of SWs by a DW, frequency doubling domain wall 
drift Ill4l3|. and magnonic spin transfer torque 
A complication that has received much less attention is 
the breakdown of the continuum approximation of mag- 
netization in materials with high anisotropy, in which the 
DW width S can be as small as a few lattice constants a 
17|. Atomic-scale DWs can display very different static 
and dynamic properties than predicted by a continuum 
model For example, discrete DW jumps that match 
the lattice periodicity have been observed by Novoselov 
et al. in thin films of yttrium- iron garnet (YIG) that sup- 
port DWs as narrow asi5/asa6/7r [19[ . Only one attempt 
to formulate the theory of SWs propagating through ul- 
tranarrow DWs is known to us [20[, predicting that a 
ferromagnetic SW can pass through a DW with little re- 
flection if its wavelength is less than twice the thickness 
of the wall which does not agree with our findings 
reported here. 

In this Letter, we address DWs in wires of insulating 
ferromagnets such a YIG that have superior magnetic 
quality factors, with damping constants that can be three 
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FIG. 1: (Golor online) Principle ground state spin configura- 
tions for a symmetric DW: its center either coincides with one 
of the lattice sites (type A) or lies between them (type B). 5 
is the domain wall width and a the lattice constant. 



orders of magnitude smaller than those in metallic mag- 
nets. Their application to novel memory concepts there- 
fore appears attractive. Charge current-induced DW mo- 
tion is ruled out, of course, but heat currents can serve 
identical purposes by virtue of the thermal spin-transfer 



torque [IJ, |21[ . Here we study SW excitations in ferro- 
magnetically coupled spin chains in the presence of a DW, 
predicting remarkable effects of the discrete lattice on 
the SW transmission. While we confirm previous results 
that in a continuum model for DWs the SW transmission 
probability is unity, spin wave reflection in a discrete spin 
model is found to be finite for abrupt domain walls. The 
resulting domain wall heat conductance should help to 
detect and manipulate DWs in insulating ferromagnetic 
wires. 

We consider a wire that is sufficiently thin such that 
the lateral degrees of freedom arc frozen out and a ID 
model is appropriate. Our model is a classical Heisenberg 
spin chain with local anisotropy and nearest-neighbor fer- 
romagnetic exchange interaction, which is appropriate for 
materials such as YIG which has a large effective spin 
(S ~ 14.3) per unit cell In the Hamiltonian [i^ 
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the first term describes the nearest neighbor interaction 
with ferromagnetic coupling J > 0, while the second 
terms is a local easy-axis anisotropy in the z-direction 
with D > 0. The local spin variable Sn is a three- 
component unit vector on the n-th lattice site. A 
transverse hard axis (shape) anisotropy causes compli- 
cated spin wave dispersions in thin films. In our one- 
dimensional system, the dipolar interaction can be im- 
portant as well, but only leads to an increased easy-axis 
anisotropy [23|, |2J] . For materials with a large unit cell 
like YIG, our approximation to integrate the magnetic 
moment distribution in a single unit cell and replace them 
with one classical spin is a good approximation at not too 
high temperatures [25| . In insulating ferromagnets, the 
heat current is associated with both phonon and magnon 
channels. In YIG the magnons, responsible for interac- 
tions with DW of interest here, contribute 2/3 of the total 
heat current at low temperatures [26 , 27 1 . We consider a 
magnetic wire of length L that connects two large heat 
reservoirs held at two temperatures T^^Tr, with constant 
difference AT = (Tl - Tr) > 0. The magnon-mediated 
heat current carried by the spin chain is given by the 
Landaucr-Biittikcr formula [28| 

Q = j dkpkVkhxJk \tk\^ [ns {i^k,TL) - ub {ujk,TR)] 

= ^ y" dujruj\tk{^^\'^ [uB iu},TL) ~ nBiuJ,TR)] . (2) 

Here pk = L/2tt is the magnon density of states, 
Vk = du/dk is the SW group velocity, ftw/c = 
2D + 2 J [1 — cos (fca)] is the SW spectrum, tk is 
the SW transmission coefficient, and nB{oJ,Tij,R) = 
Ij (^^nu/iksTi^.R) is the Bosc-Einstein distribution. 
In the absence of DWs, tk — \ and for a small tem- 
perature bias 
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where T = {Tl + Tr) /2. 

The heat current in the presence of DW is determined 
by the wall profile which is obtained by energy minimiza- 
tion with respect to the polar angle 6„ of Sn with the easy 
axis: 

sin {dn - e^-i) ~ sin (0„+i - ^n) + y sin (20„) = 0, (4) 

where we chose that boundary conditions that spins are 
oriented in the ±z directions at the ends of the wire. 
The ground state DW configuration can be computed 
iteratively [l^. To second order in the small parame- 
ter D/J<^1 this procedure yields the Walker solution 
In tan (0/2) = z/S with DW width S = a^J/ {2D) 
In a typical Walker DW, labeled A (metastable) in Fig. 
1, the spin of the central lattice site is oriented normal 
to the easy axis. However, numerical calculation leads 



to a stable type B wall (Fig. 1) which has always lower 
energies [3l ■ Clearly, this difference is immaterial in the 
continuum limit for the magnetization (micromagnetics), 
and it is only relevant for atomically sharp domain walls. 
We will discuss below that the atomic-scale magnetiza- 
tion texture is observable in the heat and spin transport 
properties. 

Proceeding from the ground state DW configuration 
we can construct, linearize, and solve the equations of 
motion dSn/dt = ~Sn x ^—ST-l/SSn^ for the spins at 
a site n in order to determine the frequencies and am- 
plitudes of the allowed SW modes. This can be done 
conveniently in a coordinate system rotated about the x- 
axis such that the transformed spins at equilibrium point 
to the new Z-axis 3l[ , where the spin vectors in the two 
coordinate systems are related as 



(5) 



The low-energy excitations correspond now to a small- 
angle precession around the Z-axis with 5f ~ 1 
The equations of motion for the SW amplitude at site n 
and frequency uj in the rotated frame read 
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where t„ = cos {On — 9n+i)- We obtain the eigenfrequen- 
cies oj and amplitudes 5'^ , numerically. 

D/ J = 0.001 leads to a broad DW profile (solid line in 
Fig. 2(a)) that cannot be distinguished from the Walker 
solution. The SW mode profile is also plotted (cir- 
cles), with the corresponding frequency given in the plot. 
For a broad DW the SW profile agrees very well with 
the analytical solution (dashed line in Fig. 2(a)) of the 
continuum model governed by the Schrodinger equa- 
tion 
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with = z/S, tp = Sx ~ iSy, q = kS with spin- wave 
vector k, and potential V {£,) ~ — 2sech^^. The poten- 
tial obtained here originates from the non-commutativity 
of t he g auge transformation with the momentum opera- 
tor |l6|. similar to electron transport in magnetic tex- 
tures [33[. However, in contrast to electrons, SWs do 
not feel a vector potential and only a scalar potential 
remains. The solutions of Eq. ([7]) are one bound SW 



D/j=6.obi ' ' ' ' ' ' 

Domain wall profile 

«oi/J=0.00649 
" Fitting 



151 




D/J=0.5 (bV 
— = — Exact domain wall profile 
- ^- ■ Walker profile 



(ico/J=1 .51031 



-150 -100 -SO 50 100 150 

z/a 

D/J=1 




- Exact domain wall profile 
■ Walker profili 
*i,i;j=2,135,27. 



(C) 0.10 




0.09 
- 0.06 
0.03 

u 

0.00 
-0.03 
-0.06 



6.4 
5.6 
4.8 

34.0; 



Bound mode 1 




Bound mode 2 




■- - -Bottom of SW band ^ 


^ * 


TopofSW band ^ ■ 




























^ Y 




V. 


, (a) 



0.4 0.62/3 O.f 

D/J 




0.2 

k (27t/a) 



FIG. 3: (Color online) (a) Frequencies of DW bound states 
and SW continuum (shaded area between dotted and dash- 
dotted lines). The solid curve represents the lowest mode, 
while the dashed curve shows the second bound state, (b) 
fc-dependence of SW transmission for different values of D/J. 



FIG. 2: (Color online) (a) The magnetization profile of a wide 
head-to- head DW (solid line). The numerical SW solution 
(circles) for frequency oj in a linear ferromagnetic chain can- 
not be distinguished from the solution (dashed line) for the 
continuum model, (b) A computed narrow DW profile (tri- 
angles) compared with the Walker model (squares). The SW 
amplitude for a selected frequency is shown by circles, (c) 
The ground state profile of an abrupt DW (triangles) com- 
pared with the Walker profile (squares). The amplitude of a 
SW for another given frequency is shown by circles. 



state and propagating SWs that arc not reflected at the 
DW but pick up a wave vector fc-dcpcndent phase shift 
T] (fc) = 2tan~^ {1/ (fc(5)} [31 that modifies the frequency 
dispersion [s^] and causes interference effects in rings [1]. 
Since the DW transmission coefficient is tk = e"'*^'^^ and 
the transmission probability \tk^ = 1, there is no ef- 
fect on DC heat transport. This conclusion is not mod- 
ified by a transverse hard-axis anisotropy D±^ {Sn)^ ^ 
since propagating SWs are elliptical but still refiection- 
less. The conclusion holds in spite of the bound magnon 
state at the DW and in contrast to conducting ferro- 
magnets, where the DW leads to an increase of electric 
resistivity even in the continuum model [sHj . In order to 
detect such a DW bound state in insulating ferromagnets, 
one has to resort to spectroscopic or other techniques [sif . 

The continuum approach breaks down for atomic scale 
DWs such as in Fig. 1. The profiles of a relatively nar- 
row DW for both exact numerical (triangles) and Walker 
solution (squares) are shown in Fig. 2(b) for a large 
anisotropy D/J = 0.5. We see that now spin waves do 
get reflected! Figure 2(c) shows the configuration of an 
abrupt DW (triangles), with D/J = 1. The Walker pro- 
file is also plotted for comparison (squares) . Interestingly, 
we now find that all SWs (circles) are totally reflected. 
The physics is explained below. 

As pointed out before, the DW supports bound SW 
states. Figure 3(a) shows the numerical results for the 



bound state frequencies as a function oi D / J below the 
continuum SW region, where the dotted line delineates 
the SW gap (fiwg = 2D) and the dash-dotted line the top 
of the SW band (fiwt = 2D + 4J). We confirm the exis- 
tence of a transition point at D/J = 2/3 that separates 
smooth and sharp DW configurations [S^] but without 
assuming \9n±i Qn\ -C 1, which is clearly not the case 
for these narrow DWs. The transition point obtained 
here differs from the value -D/ J = 5/9 calculated in Rcf . 
[3^ . Bound state 1 (solid curve), below the transition 
point Wf,i = is caused by the rotational symmetry of 
the smooth DW about the z axis. When D/J > 2/3 
the DW width is the lattice constant and bound state 1 
becomes degenerate with bound state 2 (dashed curve), 
i.e., Wfci = UJ^,2^ with frequency monotonously increasing 
with D/J. For wide DWs, i.e., D/J < 2/3, bound state 
2 merges with the continuum states. As the DW be- 
comes narrower (increasing D/J), ujb2 splits from the SW 
continuum and displays a non-monotonic dependence on 
D/J, with maximum at fiw^ ~ 0.47Z?, dropping to zero 
when D/J 2/3. 

We can now understand why the SW is reflected at 
ultra-narrow DWs. Let us consider a conflning potential 
V {£,) as 1^1 — !■ 00 for which the eigenvalue problem 
Eq. ([7]) has N bound states with q = ik„, n = 1, 2, ...,N, 
where k„ is real. Through the inverse scattering formal- 
ism, one can construct a general 2A'^-parameter formula 
for a reflectionless potential with bound states (39j , 

V{0^-'2-^lr,detA{0, (8) 

where the matrix A^n (0 = ^mn + \nK%/ + k„) , 
with Kronecker function Smn ^^^d A„ {£/) = c„e~'^"^ is 
symmetric. For A^ = 1 we recover the previous result of 
the continuum model In the discrete model, how- 

ever, N = 2. Figure 3(a) clearly shows that the bound 
state 2 vanishes into the SW continuum only for wide 
DWs [D/J < 2/3). Only in this hmit there is just 
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one bound state (mode 1) and the DW profile is per- 
fectly described by the reflectionless Walker model. To 
model the profile of a narrow DW, we expand Eq. 
to the fourth order, and obtain the DW configuration 
lntan(6l/2) = + (1/24) (a/(5)^ [6tanh (z/J) - z/(5] 

leading to the confining potential V {£,) = —2 sech^ [y (^)] 
with y (^) = ^ + (1/24) {a/Sf (6tanh^ - • According 
to Thacker et al. [39|, a reflectionless potential with two 
bound states must obey Eq. ([8]). However, the confin- 
ing potential induced by a narrow DW as derived above 
cannot be described by a four-parameter reflectionless 
potential Eq. ([S]), and indeed necessarily induces reflec- 
tions. An expansion up to higher order will not change 
this conclusion. Therefore, the SW must lose its prop- 
erty of total transmission, as confirmed by our numerical 
results. 

The fc-dependent SW transmission for different param- 
eters D/J are shown in Fig. 3(b), where the k values 
are obtained by fast Fourier transformation of the com- 
puted SW mode amplitudes jlJ]. The wave vectors k of 
interest range here from zero up to 0.87r/a correspond- 
ing to a very short wavelength (2.5a) beyond which a 
classical spin model might become unreliable [i^. In 
agreement with previous results there is no detectable 
reflection in the calculations with D/J < 0.01. How- 
ever, for a large anisotropy, the SW will be only partially 
transmitted through the DW, with transmission proba- 
bility monotonically decreasing with increasing wave vec- 
tor k. These results are very different from those ob- 
tained with a local spin-spiral approximation, i.e., a con- 
stant pitch within the DW width 4l|, which leads to 



perfect transmission for SWs with larg e fc, but hindered 
propagation at long wave lengths [20|, which is caused 
by the unphysical kinks in the magnetization texture 
of that model. We also have a message for micromag- 
netic simulations in which the continuum Hamiltonian 
n = J' {dm/dzy /2 - Dml is discretized into block 
spins J'd'^mn/dz'^ = J' (m„_i-2m„-|-m„+i) / (Az) . 
The mesh size Az is found here to be rather critical, since 
SW reflections that are artifacts in a continuum model 
occur unless 'n^pFYlD ^ 22 Az. 

Our results should be relevant for systems with close to 
atomically sharp DW such as YIG thin films. In the case 
of a narrow atomic-scale DW, the spin wave gets reflec- 
tions which leads to a modification to the heat transport 
as Q2 = Qi + AQ, where 
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We can see that the DW decreases the heat conduc- 
tance since AQ < 0. We can estimate the reduc- 
tion of the heat current due to magnon reflection at 
low temperatures as follows: magnons then exist only 
the bottom of SW band (fc = 0) at which \t{k = 0)f 
can be used to estimate the total transmission probabil- 
ity and AQ/Qi w \t{k = 0)\^ - 1. When D/J = 0.4, 



AQ/Qi ~ —35%. At room temperature fc^T is larger 
than the exchange energy J w 40 K in YIG [2^ |25| and 
magnons with high wave numbers become important for 
the heat transport. The spin wave reflection becomes sig- 
nificant when 6 < (10/7r)a, where a = 1.24 nm in YIG 

Ultranarrow DWs arc sensitive to crystal lattice pin- 
ning: the A and B-typc walls in the ID model (Fig. 1) 



have different ground state energies [18| . A domain wall 
therefore has to surmount an energy barrier when moving 
from one atomic plane to another in a "Peierls" potential 
and a threshold external force is necessary to assist the 
DW motion [l3|- Novoselov et al. [31 indeed observed 
a critical magnetic field for field-driven atomic-scale DW 
dynamics. Consequently, an intrinsic critical heat /spin 
current exists for the propagation of DWs even in per- 
fect wires. According to Ref. 15|, a temperature gra- 
dient of 1 K/nm creates a pressure corresponding to a 
field H = 5 mT for a magnetic insulator with a satu- 
ration magnetization Mg = 2 x 10^ A/m. We estimate 
the critical temperature gradient to overcome the atomic 
pinning as 0.1 K/^m for the critical magnetic field value 
= 0.7 G and = 1.6 x 10^ A/m reported in Ref. 

In conclusion, we studied the heat transport by SWs 
through magnetic DWs. When the DW becomes the 
order of a few lattice constants a second bound state 
emerges that is absent in the continuum approximation 
and causes strong SW reflection. The SW reflection leads 
to a DW heat conductance. The results developed here 
are relevant for ultrahigh-density storage devices based 
on DWs [i^l when based on insulating ferromagnets. 
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